Optimal Pointers for Joint Measurement of (t^ and cFz via 

Homodyne Detection 



Bas Janssens* and Luc Bouten^ 

*Department of Mathematics, Radboud University Nijmegen, 
Toernooiveld 1, 6525 ED, Nijmegen, The Netherlands 

'''Physical Measurement and Control 266-33, California Institute of Technology, 
1200 E. California Blvd., Pasadena, CA 91125, USA 



Abstract 

We study a model of a qubit in interaction with the electromagnetic field. By means of 
homodyne detection, the field-quadrature At + A1 is observed continuously in time. Due to the 
interaction, information about the initial state of the qubit is transferred into the field, thus 
influencing the homodyne measurement results. 'We construct random variables (pointers) on 
the probability space of homodyne measurement outcomes having distributions close to the 
initial distributions of and cr^. Using variational calculus, we flnd the pointers that are 
optimal. These optimal pointers are very close to hitting the bound imposed by Heisenberg's 
uncertainty relation on joint measurement of two non-commuting observables. 'We close the 
paper by giving the probability densities of the pointers. 



1 Introduction 



The implementation of quantum filtering and control ISl in recent experiments |13| has brought 
new interest to the field of continuous time measurement of quantum systems jll| . |12) . [S], |7], |10 |. 
123 J El- In particular, homodyne detection has played a considerable role in this development [Tn| . 
In this paper, we aim to gain insight into the transfer of information about the initial state of a 
qubit from this qubit, a two-level atom, to the homodyne photocurrcnt, which is observed in actual 
experiments. Our goal is to perform a joint measurement of two non-commuting observables in 
the initial system. In order to achieve this, we construct random variables (pointers) on the space 
of possible homodyne measurement results, having distributions close (in a sense to be defined) to 
the distributions of these observables in the initial state. 
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The problem of joint measurement of non-commuting observables has been studied by several 
authors before, see |221, ^^'^ the references therein. As a measure for the quality of 

an unbiased measurement, we use the difference between the variance of the pointer in the final 
state and the variance of the observable in the initial state, evaluated in the worst case initial 
state P^. In other words, the quality of measurement is given in terms of the worst case added 
variance. These worst case added variances for two pointers, corresponding to two non-commuting 
observables of the initial system, satisfy a Heisenberg-like relation that bounds how well their joint 
measurement can be performed 18 . 

The paper concentrates on the example of a qubit coupled to the quantized electromagnetic field. 
We study this system in the weak coupling limit jl5| . i.e. the interaction between qubit and 
field is governed by a quantum stochastic differential equation in the sense of Hudson and Par- 
thasarathy ^7]. In the electromagnetic field we perform a homodyne detection experiment. Its 
integrated photocurrent is the measurement result for measurement of the field-quadrature At+A^ 
continuously in time. Using the characteristic functions introduced by Barchielli and Lupieri jl], we 
find the probability density for these measurement results. In this density the x- and z-component 
of the Bloch vector of the initial state appear, indicating that homodyne detection is in fact a joint 
measurement of cr^ and az in the initial state. 

The goal of the paper is to construct random variables (pointers) on the probability space of 
homodyne measurement results having distributions as close as possible to those of the observa- 
bles and cr^ in the initial state of the qubit. 'As close as possible' is taken to mean that the 
pointer must give an unbiased estimate of the observable, with its worst case added variance as low 
as possible. Using an argument due to Wiseman 24 , we first show that optimal random variables 
will only depend on the endpoint of a weighted path of the integrated photocurrent. Allowed 
to restrict our attention to this smaller class of pointers, we are able to use standard variational 
calculus to obtain the optimal random variables. They do not achieve the bound imposed by the 
Heisenberg-like relation for the worst case added variances JHIj but will be off by less than 5.6%. 

The remainder of the paper is organized as follows. In Scction|21we introduce the model of the qubit 
coupled to the field in the weak coupling limit. SectionOintroduces the quality of a measurement in 
terms of the worst case added variance. This section also contains the Heisenberg-like relation for 
joint measurement. In Section 01 we calculate the characteristic function of Barchielli and Lupieri 
for the homodyne detection experiment. Section El deals with the variational calculus to find the 
optimal pointers. In Section (HI we calculate the densities of the optimal pointers and then capture 
our main results graphically. In the last section we discuss our results. 



2 The model 

We consider a two-level atom, i.e. a qubit, in interaction with the quantized electromagnetic field. 
The qubit is described by and the electromagnetic field by the symmetric Fock space T over the 
Hilbert space of quadratically integrable functions L^(R) (space of one-photon wave functions), 
i.e. 

oo 

J^:=C©0L2(R)«5=fe. 

k=\ 
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With the Fock space T we can describe superpositions of field-states with different numbers of 
photons. The joint system of qubit and field together is described by the Hilbert space ® T . 

The interaction between the qubit and the electromagnetic field is studied in the weak coupling 
limit jl4| . |15| . This means that in the interaction picture the unitary dynamics of the qubit 
and the field together is given by a quantum stochastic differential equation (QSDE) in the sense 
of Hudson and Parthasarathy ^2] 

dVt = |cr-dA* - a+dAt - ^cr+a-dt^Ut, with o-_ = ^ , cr+ = ^ ^ , Uq = I. (f) 

The operators (t_ and (t+ are the annihilator and creator on the two-level system. The field 
annihilation and creation processes are denoted At and A^ , respectively. Keep in mind that the 
evolution Ut acts nontrivially on the combined system (E) T ^ whereas (j± and At are understood 
to designate the single system-operators a±®l and I® At- Throughout the paper we will remain 
in the interaction picture. Equation |^ should be understood as a shorthand for the integral 
equation 

Ut^I+ / a-UrdAl~ / a+UrdAr~^ / a+a^Urdr, 
Jq Jo Jo 

where the integrals on the right-hand side are stochastic integrals in the sense of Hudson and 
Parthasarathy p[7|. The value of these integrals does not lie in their actual definition (on which 
we will not comment further) , but in the Ito rule satisfied by them, allowing for easy calculations. 

Theorem 2.1: (Quantum ltd rule J177. ]2Ujl) Let Xt and Yt be stochastic integrals of the form 

dXt = CtdAt + DtdA* + Etdt 
dYt = FtdAt + GtdA; + Htdt 

for some stochastically integrable processes Ct, Dt, Et, Ft, Gt and Ht (see |i7f . \2Ujl for definitions). 
Then the process XtYt satisfies the relation 

d{XtYt) = XtdYt + {dXt)Yt + dXtdYt, 

where dXtdYt should be evaluated according to the quantum ltd table: 





dA, 


dA: 


dt 


dAt 





dt 





dA; 











dt 












i.e. dXtdYt = CtGtdt. 



As a corollary we have that, for any / G C^(R), the process f{Xt) satisfies d{f{Xt)) — f'{Xt)dXt + 
^f"{Xt)(dXt)'^, where (dXt)"^ should be evaluated according to the quantum Ito table. 
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First a matter of notation. The quantum Ito rule will be used for calculating differentials of 
products of stochastic integrals. Let be stochastic integrals. Then we write 

d{ZiZ2...Zp)^ J2 M 

uC{l,...,p} 

where the sum runs over all non-empty subsets of {1, . . . ,p}. For any r/ = {n, . . . , ik}, the term [v] 
is the contribution to d{ZiZ2 ■ ■ ■ Zp) coming from differentiating only the terms with indices in the 
set {ii, . . . , ik} and preserving the order of the factors in the product. The differential d(Z\ZiZ-i)^ 
for example, contains terms of type [1], [2], [3], [12], [13], [23] and [123]. We have [2] = Zi(dZ2)Z3, 
[13] = (dZi)Z2(dZ3), [123] = (dZi)(dZ2)(dZ3), etc. 

Let us return to equation f^l- In order to illustrate how the quantum Ito rule will be used, we 
calculate the time evolution on the qubit explicitly. The algebra of qubit-observables is the algebra 
of 2 X 2-matrices, denoted M2(C). The algebra of observables in the field is given by B{T), the 
bounded operators on T. If id : M2(C) — > Af2(C) is the identity map and : B{T) — > C is the 
expectation with respect to the vacuum state $:=l®0©0©...e:F (i.e. 4>{Y) ($, 1^$)), 
then time evolution on the qubit Tt : M2(C) -> M2{C) is given by Tt{X) := id « (j){Ut*X <Si lUt). 
On the combined system, the full time evolution jt : Af2(C) (E) B{T) Af2(C) ® B{T) is given by 
jt{W) :— UfWUt- In a diagram this reads 

M2 ^* ) M2 

id(g)/ 

M2®B{T) ^' ) M2®B{T). 

In the Schrodinger picture the arrows would be reversed. A qubit-state p would be extended with 
the vacuum to p ® ij)^ then time evolved with Ut, and in the last step the partial trace over the 
field would be taken, resulting in the state poTt. 

Using the Ito rule we can derive a (matrix- valued) differential equation for Tt(X), i.e. 
dTt{X) = id ® <j){d{U;X ® lUt)) 

^id® (j){{dU*)X ® lUt + U;X ® I{dUt) + {dUl)X ® I{dUt)) 
= id (g><l){u;L{X)(g)IUt)dt 
^Tt{L{X))dt, 
where L is the Lindblad generator 

L{X) -^{a+a^X + Xa+(T^) + a+Xa^. 
In the derivation |j3Jl we used the QSDE for Uf which easily follows from Q 
dU; = U;^^a+dAt - a-dA; - ^a+a^dtj, ^ I. 

Furthermore, we used that stochastic integrals with respect to dAt and dA^ vanish with respect 
to the vacuum expectation, leaving us only with the dt terms. The differential equation ^ with 
initial condition To{X) = X is solved by Tt{X) — exp{tL){X), which is exactly the time evolution 
of a two-level system spontaneously decaying to the ground state, as it should be. Although the 
arguments above are completely standard (cf. ^7]), they do illustrate nicely and briefly some of 
the techniques used also in following sections. 
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3 Quality of information transfer 



Now suppose we do a homodyne detection experiment, enabling us to measure the observables 
At + Al in the field continuously in time If initially the qubit is in state p, then at time t the 
qubit and field together are in a state p* on M2(C) (g) given by p\W) := p (g) cpiUtWUt) = 

p{id(E)(j){U;WUt)). Since 

d(id (g) 0(C/t*/ ® {At + A;)Ut)^ = id ® (l){d{U*tI ® {At + A*t)Ut)) 

= id ® (t>{[l] + [2] + [3] + [12] + [13] + [23] + [123]) 

= lA® cj){ui{(7- +a+)(g> lUt'jdt 

— exp{tL){(T^ + (Tj^)dt = e^'^cTxdt, 

we have that regardless the initial state p of the qubit, the expectation of {At + A^) in the final 
state p* will equal the expectation of (2 — 2e~2)o-^ in the initial state p. 



3.1 Defining the quafity of information transfer 

The process at hand is thus a transfer of information about ax to a 'pointer' At + A^, which can 
be read off by means of homodyne detection. This motivates the following definition. 

Definition 3.1: (Unbiased Measurement jl8p Let X be an observable of the qubit, i.e. a 
self-adjoint element of Af2(C), and let Y be an observable of the field, i.e. a self-adjoint operator 
in (or affiliated to) B{!F). An unbiased measurement M oi X with pointer Y is by definition a 
completely positive map M : B{J^) M2(C) such that M{Y) = X. 



Needless to say, for each fixed t the map M : B{J^) M2{C) given by M {B) := iA®4){U; I®BUt) 
is a measurement of ax with pointer Y = {2 — 2e^ 2)^1(^4 A^). This means that, after the 
measurement procedure of coupling to the field in the vacuum state and allowing for interaction 
with the qubit for t time units, the distribution of the measurement results of the pointer Y has 
inherited the expectation of ax, regardless of the initial state p. However, we are more ambitious 
and would like its distribution as a whole to resemble that of ax- This motivates the following 
definition. 

Definition 3.2: (Quality |18p Let M : B{T) M2(C) be an unbiased measurement of X with 
pointer Y . Then its quality a is defined by 

^2 := sup{varpoM(l^) - Varp(X)| p e 5(M2)}, 

where S{M2) denotes the state space of M2(C) (i.e. all positive normalized linear functionals 
on M2(C)). 



This means that a^ is the variance added to the initial distribution of X by the measurement 
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procedure M for the worst case initial state p. A small calculation shows that 



VarpoAf(F) - Varp(X) = p{M{Y^) - M{Yf), 

which implies that ct^ = \\M{Y'^)-M{YY\\,w\\eie X ^ \\X\\ denotes the operator norm on M2 (C) . 
In particular this shows that is positive, as one might expect. It follows from ^21 that a equals 
zero if and only if the measurement procedure M exactly carries over the distribution of X to y. 
In short, cr is a suitable measure for how well M transfers information about X to the pointer Y . 



3.2 Calculating the quality of information transfer 



Let us return to the example at hand, i.e. M{B) = id (g) 4){UII ® BUt), with field-observable 

as a pointer for ax- Let us calculate its quality, which amounts to 
evaluating M(F^) = (2 — 2e~^)~^M {{At + )^). To this aim, we wiU first introduce some ideas 
which will be of use to us in later calculations as well. 



Definition 3.3: Let / and h be real valued functions, h twice differentiable. Let Yt be given by 
dYt = f{t){dAt + dA;), Yo = 0. For X € M2(C) we define 

Fh{X,t) -.^ id ®(j){U*X®h{Yt)Ut). 

When no confusion can arise we shall shorten Fh{X,t) to Ffi{X). 



The homodyne detection experiment has given us a measurement result (the integrated photocur- 
rent) which is just the path of measurement results for At + A^ continuously in time. Given this 
result, we post-process it by weighting the increments of the path with the function f{t) and letting 
h{y) act on the result. The following lemma will considerably shorten calculations. 



Lemma 3.4: 

dFh{X) 



Fh{L{X)) + f{t)Fh'{a+X + Xa-) + \f{tfFh>-{X) 



Proof. Using the notation below Theorem 12 . II with Zi = Ut , Z2 = I <E) h{Yt) and Z3 — Ut, we find 

dFh{X) = id® 0([1] + [2] + [3] + [12] + [13] + [23] + [123]). 

Again we will use that the vacuum expectation kills all dAt and dA^ terms. Using Theorem 12.11 
we see that after the vacuum expectation the terms [1], [3] and [13] make up Fh{L{X))dt. Since 
third powers of increments are we again have [123] = 0. From 

dh{Yt) = h\Yt)f{t){dAt +dA*t) + \h"{Yt)f{tfdt, 

we find that, after taking vacuum expectations, the terms [12] and [23] make up the second term 
f{t)Fh'{(T+X + Xa^)dt and [2] provides the last term \f{tfFh"{X)dt. □ 
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We are now well-equipped to calculate M(^{At + A^^y Choose f{t) = 1 and h{x) = . (The 
maps X i-^ X wiU be denoted 3e" hereafter.) Then M{{At + A^Y) Fj^2{I) and by Lemma rOl 

^^^^^ = 2Fx(a- -2i^3e((7_ F^2(I,Q)^0. (4) 
Applying Lemma 13.41 to Fx{(T- + cr+), we obtain 

^ ±^ = -ii^^^(a_+a+) + 2f^i(a+a_), F^(a_ + a+, 0) = 0. (5) 

Finally, Fi{a+a^) satisfies 

dFi{a+<j^) 



dt 



= -Fi{a+<j^), Fi(a+(T_,0) = (6) 



Solving 10, 10 and Q successively leads first to Fi{(j^a-) = e^^cr+fT-, then to Fx{<J^ + cr^) = 
4(e~2 — e^*)(T+(T_ and finally to Fj^2{I) = 8(e~2 — f)^cr+(7_ + t/. Consequently, the quality of 
the measurement M of ax with pointer F = (2 — 2e~2)^i(At + Aj) is given by 



1. 



1)2(7, (T_ +tl 



(2-2e-^)2 



2rT+<T_ + -IN 

(2-2e-5)2 



(2-2e-^)2 

This expression takes its minimal value 2.228 at i = 2.513, leading to a quality a = 1.493. 

The calculation above has an interesting side product. The observable M{{At + A'^)'^) depends 
linearly on ct^, indicating that in addition to information on ax, also information on az in the 
initial qubit-state ends up in the measurement outcome. Indeed, if we use as a pointer 

4(e 2 - 1)2 

then we have M{Y) = az, so that M is also a measurement of CTz with pointer Y. 

Note that the pointers Y and Y commute, i.e. measuring At + A^ via the homodyne detection 
scheme is an indirect joint measurement of ax and az- If we would also like to gain some infor- 
mation about ay, we could for example sweep the measured quadrature through [0, 27r) in time by 
measuring e"^*A( -|-e~*'^*Aj instead. In this paper however, we will restrict ourselves to continuous 
time measurement of At + A^, as additional information on ay would deteriorate the quality of ax- 
and/or Uz -measurement. The following theorem is a Heisenberg-like relation that gives a bound 
on how well joint measurements can be performed. 

Theorem 3.5: (Joint Measurement ^^ISl) Let M : B{T) M2(C) be an unbiased measure- 
ment of self-adjoint observables X E A/2(C) and X e M2(C) with self-adjoint commuting pointers 
Y and Y in (or affiliated to) B{J-), respectively. Then for their corresponding qualities a and a 
the following relation holds 

2aa>\\[X,X]\\. 
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Denote by a the quality of the cr^ measurement with the pointer Y defined in Q. Since [ax, ^z] = 
~2i<jy, the quahties a and a (corresponding to the pointers Y and Y , respectively) satisfy the 
inequality 

ad > 1. (8) 
Using similar techniques as before, that is recursively calculating Fj^4,{I) via Lemma 13.41 we find 

,2 ^ 2t-4(e-|-l)^ 
8{e-i - 1)4 {e-i - 1)2 

This expression takes its minimal value 8.836 a,t t — 2.513. This leads to a quality a — 2.973, 
which means that aa — 4.437, i.e. we are far removed from hitting the bound 1 in JSJl. However, 
there is still some room for manoeuvring by post-processing of the homodyne measurement data. 



4 The weighted path 

Let us presently return to our homodyne detection experiment. We observe Ar + A* continuously 
in time, i.e. the result of our measurement is a path uj of measurement results lOt (the photocurrent 
integrated up to time r) for At+A*. This means that we have a space il of all possible measurement 
paths and that we can identify an operator Ar + A* with the map from to M mapping a 
measurement path w g f2 to the measurement result lUt at time r. That is, we have simultaneously 
diagonalized the family of commuting operators {A^ + A* | r > 0} and viewed them as random 
variables on the spectrum fl. The spectral projectors of the operators {A^ + A*\ < t < t} 
endow with a filtration of cr-algebras S^. Furthermore, the states p'^ , defined by p'^{W) := 
p^(j>{U*WUr) provide a family of consistent measures P,- on (fi, E,-), turning it into the probability 
space (f7,Et,P). (See e.g. 0.) 

We aim to find random variables on (il, Et;P) having distributions resembling those of ax and az 
in the initial state p. In the previous section we used the random variables 

2 

YiL^) = -^^ and Y{,,) = -^^1^1—-1, 2.513 (9) 

for ax and CTz, respectively. Our next goal is to find the optimal random variables, in the sense of 
the previously defined quality. 



4.1 Restricting the class of pointers 



In our specific example, M is given by M{B) = \d® (f>{U*I ® BUr)- Note that stochastic integrals 
with respect to the annihilator Aj- acting on the vacuum vector $ are zero. Therefore, we can 
modify Ut to Zt, given by 

dZr — \ a-{dAl + dAr) — ^a+a-dr >Zr, Zq = I, 
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without affecting M Therefore, for all B e B{T), we have M{B) = id <j){U;i <Si BUr) = 
id (g) (j){Z*I (g) BZr). The solution Zt can readily be found, it is given by 

Z = ( 0^ 

' \j*e-'2r(dAr + dA*) ij- 

Note that Zt, as a matrix valued function of the measurement path, is an element of M2(C) (g) Ct, 
where Ct is the commutative von Neumann algebra generated by A^- + A*, < t < t. Moreover we 
see that Zt is not a function of all the {A^. + ^* )'s separately, it is only a function of the endpoint 
of the weighted path Yt = e~5'^((iAT- + dA*) 24 . Therefore if we define St C Ct to be the 
commutative von Neumann algebra generated by Yt , then we even have Zt G M2 (C) St ■ 

Denote by C 1-^ E[C |iSt] the unique classical conditional expectation from Ct onto St that leaves cj) 
invariant, i.e. 0(E[C \St]) = <t>{C) for all C G Ct- We can extend E[ • \St] by tensoring it with the 
identity map on the 2x2 matrices to obtain a map id(E)E[ • \St] from M2{^)®Ct onto M2{'C)® St- 
From the positivity of E[ • \St] as a map between commutative algebras, it follows that id(g)E[ • \St\ 
is completely positive. Since E[ • \St\ satisfies E[CS' \St] = E[C \St]S for all C G Ct and S G St, we 
find that id(8) E[ • \St] satisfies the module property, i.e. 

id®E[ • |5t](AiBA2) = Ai(id®E[ • \St]{B)^A2, 

for all Ai,A2 G M2{'C)®St and B G M2(C)(g)Cf. Moreover, if p is a state on M2(C), then it follows 
from the invariance of under E[ • \St] that id E[ • \St] leaves p ® 4) invariant. We conclude 
that, given p on M2(C), the map id (g) E[ • \St] from M2(C) (g) Ct onto M2(C) (g) St is the unique 
conditional expectation in the noncommutative sense of |5J that leaves p® 4) invariant. We will 
use the shorthand E^^ for id (g) E[ • \St] in the following. 

Lemma 4.1: Let C G Ct be a pointer with quality uc such that M{C) — X . Then C := E[C \St] 
is also a pointer with M (C) — X , and with quality < ac- 



Proof. Note that for all states p on M2 (C) we have 

p{M{C)) = p®4{Z-*I®CZt) ^ p®c^(z-*¥.sAl®C)Zt^ ^ p®c^{¥.s,{Z*I®CZt)) 
= p® (j){Z*I ® CZt) = p(M(C)) = p{X), 

where we used the module property and the fact that Zt is an element of M2 ® St in the third 
step and the invariance of /CKg) in the fourth step. Since this holds for all states p on M2(C), we 
conclude that M{C) = X. 

As for the variance, we note first that the conditional expectation E^^ is a completely positive 
identity preserving map. Therefore, for all self-adjoint C € Ct, we have 

EsM^C^) > (e5,(/®c))' 

by the Cauchy-Schwarz inequality for completely positive maps |23| . 
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We can now apply the same strategy as before. For all states p on Af2(C) we have 

p{M{C^)) ^ p(E)(j>(Z;i(E)C^Zt) 

= p®^(Es,{Z-*I®C^Zt)^ 

= p^cb(^z;EsAl®C^)Zt) 

ThusM(C2) > M(C2), and in particular cr2, ^ \\M{C^)~-M(C)^\\ > \\M{C'^)-M(C)^\\ = a^. □ 

This has a very useful consequence: if we are looking for pointers that record, say or az in 
an optimal fashion, then it suffices to examine only pointers in St- Instead of sifting through the 
collection of all random variables on the measurement outcomes, we are thus allowed to confine 
the scope of our search to the rather transparent collection of measurable functions of Yt . In the 
following, we will look at such pointers htiYt). We will usually drop the subscript i on /i to make 
the notation lighter. 



4.2 Distribution of Yt 

At this point we are interested in the probability distribution of the random variable Yt- Its 
characteristic function 1^ is given by 

E{k) :=Ept[exp(-z/cyt)] ^ p ® ^{u; I ® eM-ikYt)Ut) = P (^cxp^fcx) , 

so that we need only calculate -FcxpC-ifcX) (-^)- For notational convenience we will replace the 
subscript exp(— i/cse) by k in the following. Using Lemma 13.41 we find the following system of 
matrix valued differential equations: 

= -ifce 2Ffc(cr_+cr+) —Fk{I), 



dt ' ' 2 



dt 

dFk{a+<7-) fc^e-' 

= -Ffc(cr+cr_) —Fk{(J+a-), 

with initially Fk{I,0)—I, i^fc(cr_ + (T+, 0) = ct- + (7+, Ffc(CT+CT_,0) = (t+it_. 

Solving this system leads to 

We define the Fourier transform to be T{f){x) := /(fc)e*'^^dfc. Then the probability 

density of Ff with respect to the Lebesgue measure is given by -^=T{E){x) — --^p(^T(^Fk{I)) {x)j . 
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Defining p(a;) :— -^=T{Fk{I)'){x), we can write 



i.e. Yt is distributed according to a Gaussian perturbed by the matrix elements of the initial 
state p(cr- + (T+) = p{px) and /9((T+cr_) = ^p{az) + \- No information about p on ay enters the 
distribution though. To gain information about Uy we would have to change our continuous time 
measurement setup, as we discussed before. If we absorb a constant (1 — e^*)^^ in the definition 
of Yt, i.e. Yt := (1 — e"*)"' e~i{dAT + dA*), then its density becomes 

p{y) = ^ (/ + (3ty{<7- +<!+)+ PUy' - 1)^+'^-) , (10) 



where Pt ■= \/l — e 



5 Variational calculus 



In Lemma l4.ll we have shown that it suffices to consider only random variables of the form h(Yt) 
for some measurable h. In equation (|10(l . we have captured the probability distribution of Yt. All 
that remains now is to calculate the optimal h, which can be done with variational calculus. 



5.1 Optimal cTx-measurement 



We seek the function h* for which the quality cr of the pointer h*{Yt) for ct^^ -measurement is 
optimal. In other words, we need 



:= 



h^iy)p{y)dy - ( h{y)p{y)dy 
to be minimal under the restriction h{y)p{y)dy — a^. 



di 

d2 



(11) 



Now cr^ is the norm of a diagonal 2 x 2-matrix with entries di and d2 ■ Both depend smoothly on 
h, but cr^ = max{di,(i2} does not. There are three possibilities: 



I) (7^ = di in some open neighborhood of h* . To find these h* , we must minimize the smooth 
functional di and then check whether di < d2. 

II) (7^ = d2 in some open neighborhood of h* . To find these h* , we must minimize d2 and check 
whether d2 < di. 

Ill) di — d2 for h* . To find these ft.*, we must minimize d\ subject to the condition d\ = d2. 
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In principle, we need three different functionals Ai, A2 and A3 for these three distinct cases. 
However, it turns out that we can make due with the foUowing functional 

-= / h'{y)e~'^y'dy - l) + / h'{y){y^ - l)e-i^'dy) + 

71 ( / h{y)e---y'dy) + 72 ( / hiy)iy^ ~ 1)6-^^" dy) 



+ 



73 / h{y)ye-^y'dy-l). 



(12) 



The constants 71,72 and 73 are the Lagrange multipliers enforcing h{y)p(y)dy — a^- These 
are needed in all cases: Ai, A2 and A3. One can readily check that setting k = in A yields Ai, 
setting K = 1 yields A2 and considering k as a free Lagrange multiplier forces d\ = d2, so that one 
has A3 = A. 

All three cases lead to similar optimality conditions. The requirement that the optimal solution 
be stable under first order variations yields h* satisfying either 

h*i.) = ^l^±^ + C, (13) 

or 

h* (x) = dx^ + C5X + Ce (14) 
for some real constants Ci, C2, C3, C4, C5, Cg and e depending on k, 71, 72, 73- 

Suppose that h* takes the form (114(1 . The constraint h* {y)p{y)dy = a- + a+ will then force 
C4 = Ce = and C5 = P^T^ , so that h*{y) = (5^^y. The random variable we are investigating 
is simply the observed path, weighted by the function /(t) = P^^e^'^^'^, with t the final time of 
measurement. Since all the integrals we encounter are Gaussian moments, we can readily compute 
M{h*\Yt)) = !^^h*^{y)p{y)dy to be 2a+a_ + Thus 

= II (2a+a_ + P^^I) - (a_ + cj+f\\ - 1 + /3r' • 

For i — > c», this amounts to ^ V2- Already, we have improved on the naive result a = 1.493 
obtained previously. 

We proceed with the more involved case (|13|) . which will provide us with the optimal solution. 
Before we continue with the constants Ci,C2,C3 and e however, we calculate some integrals for 
later use. 



Definition 5.1: Define the error function erf(a;) and integrals /(e) and J(e) by 

2 r _ 2 e"^ /■°° e-^ 

erf(a;) —= I e ^ du, I{e) := / , , dx, J{e) := / — -2^2;. 



Lemma 5.2: 
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Proof. Since the Fourier transform of e ^l'^' is equal to \/ v^^Te' ^^'^ 



e 2 dk 



2n 1 



e du = Tr\ I — ~ erf ^ y ^ 



where, in the second step, we have used that the Fourier transform is unitary. The expression 
for J foUows from 



xe 2 



x-" + e 



dx 



dx 



\~x^ 



2x^ 



2 dx 



-1 



2e 



e - 1 



x"^ + e {x^ + eY 



x^ + s {x^ + s)^ 
dx = -V2^ + (e - l)^(e) + 2e J(e) 



□ 



The condition h* {y)p{y)dy = + = ax imphes 

V2^ 



Ci = 



^2 — C3 







A(V2^-£/(e))' 

which fixes Ci as a function of e. The next step is to express di and ^2 in terms of e: 



CI 



-,e 2 dy~l 



CI 



di 



277 J -00 iy^ + 



27r 7_oo (y' + £)' 



e 2 (fiy 



= (l{e)-eJ{e))-l, 

^(%/2^-(l + 2e)/(e)+e(l + e)J(e)) +^ 



First, we use Lemma 15.21 to express the above in terms of elementary functions and the error 
function. Then, using Maple, we find that e 1— > max{c?i, ^2} has a unique minimum at e = 0.605, 
for which di = d2 ~ 0.470. This leads to a Ci that equals 2.359, and to a quality of 



a — A/max{di, ^2} — 0.685 . 



5.2 Optimal a^-measurement 

For optimal a^-measurement, we can run the same program. We search for the function h which 
optimizes the quality (t, under the restriction that h{Yt) be a pointer for cr^ -measurement. That is, 
we search for a function h minimizing the functional of equation , but now under the restriction 
/oo ^{y)p{y)'^y — Again there are three cases of interest, di — d2, di > c?2 and o?2 > ^i, which 
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we can treat simultaneously by introducing, analogous to equation (|12|) . the functional 



A(/i,«;,7i,72,73) := / h^y)e-^y"dy-l)+t,(^ / h^{y){y^ - l)e-^y' dy) + 

71 J^^h{y)e~^y'dy+l) +72(^ j^^h{y){y^ - l)e-'^y' dy - 2) + 

. /-oo 

73 ( / h{y)ye^^y^dyj. 



Indeed, 71,72 and 73 are the Lagrange multipliers enforcing the restriction h{y)p{y)dy — dz- 
Again, the functional of equation depends non-differentiably on h when di = d2- We then 
have to search for the optimum among the points of non-differentiability, in which case n is the 
Lagrange multiplier confining us to these points. If di > ^2 then k = 1 and if ^2 > di then k = 0. 
Summarizing, wherever A takes its minimal value, optimality implies k, 71, 72, 73) — for 

some K, 71,72 and 73. Performing the functional derivative yields either 

or 

h{x) = Dix"^ + D5X + Dq (16) 
for some (time-dependent) constants Di, D2, D^, D4, D^, Dq and S depending on re, 71, 72 and 73. 

Again, we begin with the least demanding case (|16|) . resulting from k 0. The condition 
h{y)p{y)dy = implies Dc, = 0, D4 — fi^"^ and ~ —1 — j3f . For t 00, this leads 

to°° 

= \\M{h\Yt)) - M{h{Yt)f\\ = ||(4r7+a_ +3/) -/|| =6, 

so that a \/&. 

This improves the result a — 2.973 obtained previously, but once again the ultimate bound will be 
reached in the more arduous case (|15|l . There, the condition h{y)p{y)dy ~ Oz implies 

This leads to expressions for di and ^2 as a function of 5. Using Lemma 15.21 and Maple once 
more, we find that the function S t—f max{(ii,d2} has a unique minimum at S ^ 2.701, for which 
di = ^2 = 2.373. This leads to a quality of 



a = •\/ max{(ii, ^2} = 1.540, 
attained for Di = 0, D2 = —21.649 and = 5.391. For the joint measurement this leads to 

o-CT = 1.056. 

Although we did not achieve the bound of 1 provided by Theorem 13. 51 we have come as close as 
the measurement setup allows. We conclude that, using the setup investigated in this article, no 
simultaneous measurement of and will be able to approach the quantum bound by more 
than 5.6 %. Furthermore, we have identified the unique pointers for this optimal measurement in 
equations ((T^ and ((T5|) . 
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6 Distribution of pointer variables 



We have designed pointers h* (Yt) and h{Yt) in such a way that their distributions in the final state 
best resemble the distributions of and cr^ in the initial state. We will now calculate and plot 
these final densities. 



6.1 Calculation of h*- and /i-densities 



Let p be the initial state of the qubit and let it be parameterized by its Bloch vector {P^, Py,Pz). 
By equation H10(l . the density q{y) of Yt is given by 

qiy) = p{p{y)) = ^ (i + PtyP. + PUy' - 1)^^) • (17) 



We are interested in the the distributions r{x) and s{x) of h*{Yt) and h(Yt) respectively. Let us 
start with h* . From equation we first calculate the points y where h*(y) = x for some fixed 
value of X. 

Ci ± JC} - 4a;2£ 

y± = -7. 

2x 

By the Frobenius-Peron equation (see e.g. ^|), r(x) is given by 

}t\h*\y±)\' 

which leads immediately to 

where it is understood that r{x) ^ only for x G [— ^fe]' 

We run a similar analysis for s{x). The points y in which h[y) = x are given by 



y± = ±1 

This leads to 



i{x-D^)5-D2 



Ds-x 



six 



(2/|+<5)2(l + /?,y±P.+/3?(y|-l) 



with s{x) 7^ only for x G [Z^a + D3]. We proceed with a graphical illustration of the results 
obtained so far. 



w± 
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6.2 Plots of cr^-measurement 

According to formula El the distribution of the endpoint of the weighted path depends on the 
input qubit-state. For instance, the negative CT^j-eigenstate, the tracial state and the positive ct^- 
eigenstate lead to the distributions below: 






Figure 1: Probability density of Figure 2: Probability density of Figure 3: Probability density of 
the endpoint of the weighted path the endpoint of the weighted path the endpoint of the weighted path 
for input | <— ). for input tr. for input | — >). 



In order to estimate ax , we use the pointer of Ux given by 

Cix 



h*{x) 



(20) 



with e = 0.605 and Ci — 2.359. It is illustrated in figure 4 
to the right. 

In formula ((TH|l . we have calculated the probability distribu- 
tions of this pointer under the distributions of the endpoint 
of the weighted path illustrated in figures 1, 2 and 3. They 
are illustrated in figures 5, 6 and 7 below. 





Figure 4: Pointer for cfx 



5-^ 




4— 




3— 




2— 

X. 1 — 




1 1 1 1 1 1 1 1 1 1 1 n 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



I I I I I I I I I I I rtl I 
-1.5 -1.0 -0.5 0. 



I I I I I I I I I I I I I 
.0 0.5 1.0 1.5 



I I I I I I I I I I I I I 
-1.5 -1.0 -0.5 0. 



I I I I I I I I I I I I I I 
,0 0.5 1.0 1.5 



Figure 5: Probability density of Figure 6: Probability density of Figure 7: Probability density of 
the (Tj^-pointer for input | <— ). the cra;-pointer for input tr. the (Ti,-pointer for input | 
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6.3 Plots of (j^-measurement 



We repeat this for the cr^-pointer. By formula [T7I the positive (Jz-eigenstate, the tracial state and 
the negative CTz-eigenstate lead to the distributions of the endpoint of the weighted path shown 
below: 





1 2 3 4 5 




Figure 8: Probability density of Figure 9: Probability density of Figure 10: Probability density of 
the endpoint of the weighted path the endpoint of the weighted path the endpoint of the weighted path 
for input | t ). for input tr. for input \ I). 



In order to estimate (Tz , we use the pointer of cr^ illustrated 
here to the right. It is given by 



h{x) 



Do 



(21) 



with S = 2.701, D2 = -21.649 and D3 = 5.391. 



Figure 11: Pointer for (Tz. 

From formula H19|l . we read off the probability distributions of this pointer under the distributions 
of the endpoint of the weighted path illustrated in figures 8, 9 and 10. They are shown below: 





Figure 12: Probability density of Figure 13: Probability density of Figure 14: Probability density of 
the (Tz -pointer for input | t )• the -pointer for input tr. the -pointer for input | J, ). 
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7 Discussion 



In this paper, we have investigated homodyne detection of spontaneous decay of a two-level atom 
into the electromagnetic field. We have seen how the photocurrent, besides carrying information on 
ax (which is immediate from the innovations term in the filtering equation), also carries information 
on Gz- Homodyne detection can thus be viewed as a joint measurement of the non-commuting 
observables Gx and tr^ in the initial state of the qubit, and we have identified the optimal pointers 
for this procedure in equations and (|^ . 

One particular feature of the pointers we constructed might seem counterintuitive at first: they 
yield values outside [— 1 , 1] with nonzero probability. This is a direct result of our requirement that 
the measurement be unbiased. Suppose, for example, that the input state is | so that has 
value 1 . Since the photocurrent carries information on Ux as well, its information on az is certainly 
flawed, and will yield estimates tr^ < 1 at least some of the time. Unbiasedness then implies that 
also estimates cr^ > 1 must occur. 

On the other hand, an unbiased measurement will yield on average the 'true' value of for any 
possible input state. (Not just for the 3 possibilities sketched on page llTI ) In repeated experiments, 
optimality of our pointers guarantees fast convergence to these averages. 

Theorem 13.51 provides a theoretical bound for the quality of joint measurement of Gx and cr^. No 
conceivable measurement procedure can ever achieve a a < 1. It is now clear that this bound 
cannot be met by way of homodyne detection: a small part of the information extracted from 
the atom is simply lost in this particular procedure. Constructing the optimal pointers on the 
photocurrent does yield aa = 1.056 however, a figure much closer to the bound than the 4.437 
provided by the naive choice of ®. 
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